Approximate action-angle variables for the figure-eight and other periodic three-body 

orbits 
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We use the maximally permutation symmetric set of three-body coordinates, that consist of the 
"hyper-radius" R — \J p 1 + A 2 , the "rescaled area of the triangle" -£0i\p x A|) and the (braiding) 
hyper-angle <j> — arctan ^ \i^pi ) i to analyze the "figure-eight" choreographic three-body motion 

discovered by Moore 0] in the Newtonian three-body problem. Here p, A are the two Jacobi relative 
coordinate vectors. We show that the periodicity of this motion is closely related to the braiding 
hyper-angle (j>. We construct an approximate integral of motion G that together with the hyper- 
angle (f> forms the action-angle pair of variables for this problem and show that it is the underlying 
cause of figure-eight motion's stability. We construct figure-eight orbits in two other attractive 
permutation-symmetric three-body potentials. We compare the figure-eight orbits in these three 
potentials and discuss their generic features, as well as their differences. We apply these variables 
to two new periodic, but non-choreographic orbits: One has a continuously rising <f) in time t, just 
like the figure-eight motion, but with a different, more complex periodicity, whereas the other one 
has an oscillating (f>(t) temporal behavior. 

PACS numbers: 45, 45.50.-j, 45.50.Jf, 5.45.-a 



I. INTRODUCTION 

The three-body problem is one of the oldest and most 
challenging in classical mechanics pj • Until recently only 
a few periodic three-body solutions were known [l| in 
Newton's gravitational interaction potential. A new pe- 
riodic, "figure eight", trajectory was found in 1993 by 
Moore Q in the case of three equal masses and gravita- 
tional — 1/r potential, using numerical methods. Its ex- 
istence and stability were later proven formally by way of 
variational arguments [3j, but no closed (analytic) form 
of this solution has been shown as yet. Moreover, the 
figure-eight solution has also been found in the general- 
relativistic three-body dynamics [9], and its bifurcations 
have been studied as a function of the mass asymmetry 
Q. Proofs of existence, as well as some properties of 
figure-eight orbits in pairwise sums of — l/r a two-bod; 
potentials with 1, have been studied in Refs. 
[7|. Any new solution, and/or insight into the existing 
ones should be of intrinsic interest. 

Of course, the figure-eight orbit is highly symmetric, 
but it is not immediately clear what is the underlying 



'Electronic address: 
t Electronic address: 



suki@ipb.ac.rs 



dmitrasin@ipb.ac.rs 



dynamical reason for this symmetry. It is an empiri- 
cal fact, however, that all known figure-eight orbits exist 
only in (three-body) permutation symmetric potentials. 
Indeed, it is known that the figure-eight orbit bifurcates 
into new, less symmetric orbits as one changes the mass 
ratio(s) of the three particles, and thus breaks the per- 
mutation symmetry, see Refs. Q, We explore this 
connection between the permutation symmetry and the 
figure-eight orbit and make it more explicit. In the pro- 
cess we have found new solutions with lesser symmetry, 
much like those in Ref. [18} , and obtained new insights 
into the role of permutation symmetry in the classical 
three-body problem. 

In this paper we report our studies of figure-eight or- 
bits in three kinds of three-body potentials: 1) the New- 
tonian gravity, i.e., the pairwise sum of —1/r two-body 
potentials; 2) the pairwise sum of linearly rising r two- 
body potentials (a.k.a. the A string potential); 3) the 
Y-junction string potential [Hj],[l6| that contains both 
a genuine three-body part, as well as two-body contri- 
butions (this is the first time that the figure-eight has 
been found in these string potentials, to our knowledge). 
These three potentials share two common features, viz. 
they are attractive and symmetric under permutations of 
any two, or three particles (2fjj . 

A set of variables makes this permutation symme- 
try manifest and we use them to plot the motion of 
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a numerically calculated figure-eight orbit. As there 
are three independent three-body variables, and there 
can be at most two independent permutation-symmetric 
three-body variables [2l|, the third variable cannot be 
permutation-symmetric. In other words the third vari- 
able must change under permutations. Moreover, it must 
be a continuous variable and not be restricted only to 
a discrete set of points, as is natural for permutations. 
Thus it must provide a smooth interpolation between 
(discrete) permutations. We identify here the third inde- 
pendent variable as 4> — arctan ( ^fj^i J and show that 

it grows/descends (almost) linearly with the time t spent 
on the figure-eight trajectory and reaches ±27r after one 
period T. Thus, <fi is, for most practical purposes, inter- 
changeable with the time variable t on the figure-eight 
orbit. The hyper-angle </> is the continuous braiding vari- 
able that interpolates smoothly between permutations 
and thus plays a fundamental role in the braiding sym- 
metry of the figure-eight orbits 0, ■ 

Then we construct the hyper-angular momentum G3 = 
i (p p • A — pa • p) conjugate to <f>, the two forming an 
(approximate) pair of action-angle variables for this pe- 
riodic motion. Here we calculate numerically and plot 
the temporal variation of </>, as well as that of the hyper- 
angular momentum Gs(i), the hyper-radius R(t) and 
r(t). We show that the hyper-radius R(t) oscillates about 
its average value R with the same angular frequency 
(3<p) and phase, as the new ("reduced area") variable 
r(t). Thus, we show that <fi(t) is, for most practical pur- 
poses, interchangeable with the time variable i, in agree- 
ment with the tacit assumption(s) made in Refs. @, 
[H , though the degree of linearity of this relationship de- 
pends on the precise functional form of the three-body 
potential, see Sect. IIIIBI 

As stated above, </> is not exactly proportional to time 
t, but contains some non-linearities that depend on the 
specifics of the three-body potential; consequently the 
hyper-angular momentum G3 is not an exact constant 
of this motion, but oscillates about the average value 
G3, with the same basic frequency 30. Thus, the time- 
averaged hyper-angular momentum G3 is the action vari- 
able conjugate to the linearized hyper-angle 4> ■ 

We use these insights to characterize two new planar 
periodic, but not choreographic three-body motions with 
vanishing total angular momentum. One of these orbits 
corresponds to a modification of the figure-eight orbit 
with (j)(t) that also grows more or less linearly in time, 
but has a more complicated periodicity pattern defined 
by the zeros of the area of the triangle formed by the 
three particles (also known as "eclipses" , "conjunctions" 
or "syzygies"). Another new orbit has (f>(t) that grows 
in time up to a point, then stops and "swings back". 
We show that this motion, and the other two, can be 
understood in view of the analogy between the three- 
body hyper-angular ("shape space") Hamiltonian on one 
hand and a variable-length pendulum in an azimuthally 
periodic in-homogeneous gravitational field, on the other. 



This paper is divided into five parts: after the Intro- 
duction in Sect. [IT] we introduce a complete (maximal) 
set of permutation symmetric three-body variables and 
illustrate them with two examples: 1) the curves in the 
"shape space" of triangles depicting those triangles with 
one of its three angles equal to a particular value in the 
range (f , 7r); 2) the contour plots of the Newtonian grav- 
ity, the Y-junction string and the A-string potentials. In 
Sect. ITTTJwe show the time dependence of the figure-eight 
motion in Newton's gravity and the Y-string potentials. 
In Sect. IIVI we show and discuss the new solutions. Fi- 
nally in Sect. [V] we summarize and draw conclusions. 



II. PERMUTATION SYMMETRIC 
THREE-BODY COORDINATES 

As the static three-body potential depends on three in- 
dependent scalar variables, e.g. the pairwise relative dis- 
tances/separations, the choice of appropriate (relative) 
variables is a crucial one. A number of three-body rela- 
tive variables have been devised, starting with those in- 
troduced by C.G. Jacobi in the 19th century [l(|, and 
extending to the so-called hyper-spherical coordinates in- 
troduced in the 1960's [nf, QJ], 0. These variables 
were introduced in attempts at solving certain quantum 
mechanical three-body problems, that demand special at- 
tention to be paid to the permutation symmetry. Nev- 
ertheless, only one, Ref. of these sets is manifestly 
permutation symmetric and yet it has not been widely 
used. 

Here we use the manifestly permutation symmet- 
ric three-body variables, apparently first introduced by 
Hopf: the hyper-radius R, the "scale-invariant area" of 
the triangle \J\ — r 2 = 2R~ 2 \p x A|, where, and find 

as the hyper-angle </> = arctan ( ^fj^i ^ , that is conju- 
gate to the generalized hyper-angular momentum G3 — 
\ (p p • A — pa p) ■ One may relate these to the hyper- 
spherical variables x = ^jgr- and z = A ~f that have 
the circle with unit radius as their natural domain. Then 
the area of the triangle -^|p X A| and the hyper-radius 
R are related to the the new variables r, (f> as follows 

(,-.+,■.)_!_ (3«g*iy (1) 

<t> = tan" 1 ^yj . (2) 

The hyper-angle 4> is zero at the (x = 0, Z = 1) point 
("12 o'clock") and increases as one moves clockwise. 



A. The shape space of triangles 

The natural domain of the permutation symmetric 
variables is a circle with unit radius, see Fig. [TJ The 
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FIG. 1: (Color online) The curves/lines in shape space of tri- 
angles depicting triangles with one fixed angle: a) the outer 
blue (dashed) line, for the fixed angle equal to 109.5°, b) 
blue (full straight) lines for the fixed angle equal to ^, as 

functions of z = z = cos 2\ (ordinate = vertical axis) and 
x — x\/\ — z 2 = cos 6 sin 2\ (abscissa = horizontal axis). 
The domain of these variables is a magenta (dark gray) cir- 
cle of radius unity. The two straight red (dashed) lines at 
angles of and the vertical axis are the symmetry axes, 

i.e. S2 subgroups of the S3 permutation group, and of the 
"constant angle curves" in shape space, as well. The three 
collinear configurations in which one pair of particles has van- 
ishing separation are denoted by big solid circles, and the 
three collinear configurations in which one particle has equal 
separations from the other two are denoted by small solid 
circles. 



points on the unit circle correspond to collinear configu- 
rations ("triangles" with zero area). 

The two straight lines at angles of ±^f , together with 
the vertical axis are the three (reflection) symmetry axes; 
these reflections correspond to the three "two-body per- 
mutations" /transpositions in the S3 permutation group. 
The two cyclic permutations of the S3 permutation group 
correspond to the rotations through ±^p 

The six points where the symmetry axes cross the big 
circle in Fig. Q] correspond to either a) three collinear 
configurations ("shapes") in which one pair of particles 
has vanishing separation (big solid circles), i.e. "sits on 
top of each other", or b) three collinear configurations 
("shapes") in which one particle has equal separation 
from the other two, i.e. "sits in the middle between 
the other two" (small solid circles). The center of the 
circle corresponds to the equilateral triangle configura- 
tion ( "shape" ) , which turns into a point when the hyper- 
radius R — » 0. 



B. Newton's, A and Y-string potentials 

In Figs. I2l3l4l we show three attractive three-body po- 
tentials that are either pairwise sums of two-body terms 



FIG. 2: (Color online) The equipotential contours for the cen- 
tral Y-string potential, and the boundary between the central 
Y-string and two-string potentials as functions of z = z = 
cos2x (vertical axis), and x = x\J\ — z 2 = cos#sin2x (hor- 
izontal axis). The blue (dashed) curve denotes the boundary 
between the two-body and the three-body components of this 
potential, see Ref. [15|]. The rotation symmetry about the 
axis pointing out of the plane of the figure should be visible 
to the naked eye. 



viz. Newton's 



x ^ 1 

VNowton --5 2-, | x . -x-l' 



and the A-string 



Va = Q"a ^ |x; - Xj 



(3) 



(4) 
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or contain such a two-body component in a limited part 
of the configuration (shape) space, such as the Y-string 

3 3 
Vy = cry ruin |xj — x| = (Ty |xi — xt I, (5) 

i=r i-i 

where the minimum of the sum occurs at the Torricelli 
point x = xt, see Ref. [l6j]. Note that the Y-string 
potential has perfectly concentric contour lines within a 
pear-shaped region of shape space delineated by the blue 
dashed line in Fig. [2] As shown in Ref. 1J] , that "hyper- 
rotational" symmetry leads to a new constant of motion 
in this part of shape space. A clear discrimination of the 
Y-string from the A-string three-quark potentials had 
been a problem in lattice QCD until Ref. showed 
that the two kinds of potentials have essentially different 
hyper-angular dependencies. The separation of one kind 
of three-body potential from another is facilitated by the 
use of the new variables r and <j>. Then the Y-string com- 
ponent is manifested through the sole dependence on r, 
whereas the A-string is manifested through the depen- 
dence of the potential on the hyper-angle <fi, within the 
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FIG. 3: (Color online) Contour plot of the A-string potential 
as a function of 2 = z = cos 2\ (vertical axis) and x = 
— cos 9 sin 2\ (horizontal axis) for any fixed value 
of the hyper-radius R. The rest of the legenda is as in Figs. 
[T]and[2] The center of the circle is the point with the highest 
value of the potential. 
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FIG. 4: (Color online) Contour plot of the logarithm of the 
sum of Newton's two-body potentials as a function of z = 
z — cos2x (vertical axis) and x = xyl — z 2 = cos#sin2x 
(horizontal axis) at any fixed value of the hyper-radius 7?. The 
rest of the legenda is as in Figs. [T]and[2] As one approaches 
the two-body collision points (three large solid points on the 
big circle), the equipotential contour lines become more and 
more dense and parallel, finally reaching infinite density at 
these points, due to the singularities/poles present. 

confines of the "central potential" boundary in terms of 
"old" variables ( x ,9)) 0. 

Note that all three potentials in Figs. I2I3I4I have es- 
sentially (topologically) the same form in the two hyper- 
angular ("shape space") variables: 

V(r,4>) = V{r) + 8V(r) cos(30) + . . . . (6) 

This is a consequence of their permutation symmetry. 
Any attractive permutation symmetric potential has its 



highest value at the center of the circle (r — 0)and it 
decreases monotonically as one moves radially towards 
the r — 1 circle. Moreover, a permutation symmetric po- 
tential is circularly symmetric at the center (8V(0) = 0) 
and is increasingly broken by a periodic (f> angular ( "two- 
body") component 6V(r) cos(3</>) as one moves radially 
towards the r = 1 circle. 

As a consequence of this "topological equivalence" 
these potentials lead to certain kinds of orbits, such as 
the "figure eight" one, that are essentially identical. The 
hyper-radial part of the potential does not seem to be 
very important, so long as it is attractive, because the 
stability of the orbit is ensured by the approximate (dy- 
namical) 0(2) symmetry of these potentials. The details 
of these potential differ, of course, and therefore lead to 
different detailed properties of the amplitude and phase 
variations, but not so for the qualitative properties of the 
motion. Indeed, if the potential does not contain the peri- 
odic 0-dependent "two-body" component near the "outer 
edge" of the shape space circle (or near the equator of the 
shape space hemisphere), then there is no "figure-eight" 
orbit in that potential. 

C. Approximate dynamical 0(2) symmetry 

The sum of Newton's or A-string two-body potentials 
is approximately symmetric under infinitesimal rotations 
in the shape space, at least in the central (r ~ 0) part of 
the "shape space" , as can be seen in Figs. G2 S whereas 
the Y-string is exactly symmetric in the same region, see 
Fig.H 

Of course, the two-body potentials are exactly in- 
variant under the finite ("kinematic") rotations through 
<f> = that correspond to cyclic permutations, as well 
as under reflections about the three symmetry axes, that 
correspond to binary/two-body permutations ("transpo- 
sitions"). 

Independence of the potential on the variable <f> is 
equivalent to its invariance under (infinitesimal) "kine- 
matic rotation" 0(2) transformations 

Sx = 2ez (7) 
5z = -lex, (8) 

or, in terms of the original Jacobi variables, 

Sp = eX (9) 
SX = sp. (10) 

in the six-dimensional hyper-space. This invariance leads 
to the new integral of motion G3 = | (Pp ■ A — pa • p), 
associated with the dynamical symmetry (Lie) group 
0(2) that is a subgroup of the (full hyper-spherical) 0(6) 
Lie group. 

This 0(2) symmetry transformation is an infinitesimal 
version of the so-called "kinematic rotations", see Ref. 
that operate in two, ordinarily different planes at 
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FIG. 5: (Color online) Real space trajectories of the figure- 
eight and one new solution that passes through the figure- 
eight initial configuration for three different potentials. The 
legenda are explicitly shown in the upper left-hand corner 
of the figure: 1) Newtonian potential figure-eight: red (dark 
gray solid) curve; 2) Y-string potential figure-eight: (dark 
gray long dashed) curve; 3) Y-string potential new solution: 
blue (gray medium-length dashed) curve; 4) A-string poten- 
tial figure-eight: magenta (light gray short dashed) curve. 
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FIG. 6: (Color online) The time dependence of the hyper- 
angular radius r red (solid), and the hyper-angles a = sin - r 
blue (long dashed) and cj> gray (short dashed) of the figure- 
eight solution in Newton's potential. The legenda are explic- 
itly shown below the figure. 

particle in the middle of the other two. The initial veloc- 
ities are such that the total angular momentum vanishes. 



the same time(22j: 1) in the plane of Jacobi vectors p, A, 
and 2) in the plane of Jacobi momenta p p — p\ (these 
two planes need not coincide in general). It is only in 
the special case of planar motions that these two planes 
coincide, and it is only in the (even more special) case of 
vanishing (total) angular momentum that the new con- 
stant of motion has presently discernible consequences. 

In the case of the sum of two-body potentials, such 
as the Newtonian gravity, or the A string potential, this 
generalized hyper-angular momentum G3 is not an exact 
integral of motion, but an approximate one. The precise 
consequences of such an approximate symmetry depend 
on the initial conditions of motion, as we shall see below. 



III. FIGURE-EIGHT MOTION 

A periodic "figure eight" orbit, Fig. [SJ with vanishing 
total angular momentum (L—0) has been found by Moore 
in the case of equal masses and gravitational potential. 
It was shown in Ref. [5] that the hyper-radius R is close 
to being constant along the figure-eight trajectory: it 
only makes small-amplitude oscillations in lockstep, i.e. 
with the same frequency and locked in phase, with the 
area of the triangle and the hyper-angle <p, see below. 



A. Time dependence of the hyper-angular motion 

We take (rinit., 0init.) = (1, \' K ) a s the initial condi- 
tion, which is one of three identical configurations, up to 
permutations. This is a collinear configuration with one 



1. Newtonian gravity 

In Fig. [6] we see that both hyper-angular variables 
(r, </>) oscillate with the same frequency and locked in 
phase along the figure-eight trajectory. 

2. The Y- and A string potentials 

A similar situation is present in the other two poten- 
tials: the hyper-radius R{t) is almost constant along this 
trajectory: it makes small-amplitude oscillations in phase 
with the area of the triangle, Fig. [7] Similarly for the A 
string, see Fig. [5] 

B. Hyper-angular (f> dependence 

One can see in Fig. [9] that the periodicity of the figure- 
eight motion is determined by the braiding angle <j). Here 
one can also see that the actual path in the shape space, 
Fig. [SI taken by the Newtonian three-body system is 
remarkably close to the Newtonian iso-potential lines in 
Fig. If this were true, then the hyper-radius would 
be constant along the figure-eight orbit, but it is not: 
R{t) and r(t) are periodic functions, with the same basic 
frequency of 3<t>(t) and locked in phase, i.e. they oscillate 
about their average values as follows 

4> = + <tysin(30) + ... 
r(cj>) = r + Sr sin(30) + . . . 
R((f>) = R + 5Rsm(3(t>) + .... (11) 
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FIG. 7: (Color online) The time dependence of the hyper- 
angular radius r red (solid), and the hyper-angles a = sin - r 
blue (long dashed) and <f) gray (short dashed) of the figure- 
eight solution in the Y-string potential. The legenda are ex- 
plicitly shown below the figure. 
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FIG. 8: (Color online) The time dependence of the hyper- 
angular radius r red (solid), and the hyper-angles a = sin -1 r 
blue (long dashed) and cf> gray (short dashed) of the figure- 
eight solution in the A-string potential. The legenda are ex- 
plicitly shown below the figure. 



This phase- and frequency locking provide an important 
constraint that effectively reduces the number of inde- 
pendent degrees of freedom to two. In other words (f> 
is the (approximate) cyclic, or "ignorable" variable of 
the figure-eight periodic motion that may be integrated 
out /ignored/. The conjugate action variable G3 is the 
associated (approximate) integral of motion. 



C. The action variable conjugate to <f> 

In order to find the appropriate action variable we look 
at the ( "kinematic" ) hyper-angular momentum G3 that 




FIG. 9: (Color online) Trajectories of the figure-eight and one 
new solution in three different potentials in terms of 2 = z — 
cos 2;^ (vertical axis) and x = x\/l — z 2 = cos 8 sin 2\ (hori- 
zontal axis). The legenda are explicitly shown in the middle 
of the figure: I) Newtonian potential figure-eight: red (dark 
gray solid) curve; 2) Y-string potential figure-eight: (dark 
gray long dashed) curve; 3) Y-string potential new solution: 
blue (gray medium-length dashed) curve; 4) A-string poten- 
tial figure-eight: magenta (light gray short dashed) curve; 5) 
unit circle: magenta (dark gray solid). 



reads 

G 3 = ^(i? r ) 2 0= jiRsinafj) (12) 

(with vanishing angular momentum L = 0) as a function 
of permutation symmetric variables R,r = sin a and <f>, 
and oscillates as a periodic function of the hyper-angle 
3(f). Hence it follows that G3 and <j) must be (almost) 
constant in orbits with L = 0, R ~ const and r ~ 1. In 
other words, the angle <fr grows (or decreases, depending 
on the orientation of the motion) almost linearly in time, 
which is confirmed by our numerical results. 
The time/hyper- angle average G3 




is a non- vanishing constant on the figure-eight orbit, fur- 
nishing the (approximate) action variable that goes to- 
gether with the (linearized) hyper-angle <f> for this peri- 
odic motion. The approximate constancy of G3 ~ G3 ^ 
is the cause of dynamical stability of the figure-eight or- 
bit: The vanishing angular momentum (L — 0) three- 
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FIG. 10: (Color online) The time dependence of the hyper- 
radius R red (solid) , and the hyper-angular momentum G gray 
(dashed) of the figure-eight solution in Newton's potential. 
The legenda are explicitly shown below the figure. 
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FIG. 11: (Color online) The time dependence of the hyper- 
radius R red (solid), and one third of the hyper-angular mo- 
mentum G/3 gray (dashed) of the figure-eight solution in the 
Y-string potential. The legenda are explicitly shown below 
the figure. 



body relative motion kinetic energy 



„ III, I -o 

Tidn = - I R 2 



R 2 



1-r 



r <f> 



(j) sina^ 



(14) 



has the form of the single-particle kinetic energy in 
polar coordinates, albeit with polar radius re- 
duced by half in the "hyper-angular kinetic energy" 



sin cxj 



This means that another 



(hyper-) angular momentum-like three vector G is con- 
served when the potential does not depend on the two 
angles (a, 0). 

In the case when the potential depends on a, but does 
not depend on 0, or has only small variations with 0, then 
the "azimuthal hyper-angular momentum" G3 = 9 ^ in is 
approximately constant: 



dV(r, (j)) 
d4> 



38V(r) sin(3<£) 



— 2G 

Then G3 ~ G3 7^ provides a repulsive term 

in the effective hyper-radial potential V c g(R) — + 
V3-body(-R) that prevents the system from collapsing to 
a point, just as the (ordinary) angular momentum L ^ 
does in the two-body problem. 



1. Newtonian gravity 

The temporal variation of the hyper-radius R(t) and 
the hyper-angular momentum G(t) in the Newtonian 
gravitational potential are shown in Fig. [TOT 



2. The Y- and A string potentials 

The temporal variation of the hyper-angular momen- 
tum G(t) and R(t), the former reduced by factor three, so 
as to emphasize the small variation of R in the Y-string 
potential are shown in Fig. [TT] Note the cut off peaks of 
the sines ("flat tops") of the hyper-angular momentum 
G/3 due to the exact dynamical 0(2) symmetry in that 
part of the configuration space. 



IV. NEW Y-STRING PERIODIC ORBITS 

Note, moreover, that the permutation symmetric 
three-body potential V3_body(sina, (j>) in the region of 
the figure-eight orbit (i.e. on the outer fringes of the 
(r = sin a, (f>) circle) is attractive as a function of r — 
sin a, with a minimum at the unit circle (r = 1), and a 
strictly periodic function of the (triple) hyper-angle 3(f>. 

Then it should be no surprise that the figure-eight 
motion of the three-body system, with its almost con- 
stant hyper-radius R, has many similarities with that 
of the spherical pendulum in an inhomogeneous (az- 
imuthally periodic) gravitational field: figure-eight orbit 
corresponds to rotations, but there are other qualitatively 
different kinds of motions that we shall display and briefly 
discuss in this section. 

There is a small, yet pronounced non-linearity in 
the figure-eight motion's (^'s temporal dependence, par- 
ticularly near the </> — 0,±^ points. These three 
points/lines in the (r, </>) circle, see Fig. [9l correspond to 
the configurations of closest two-body approach in real 
space. Of course, the figure-eight orbit does not touch 
the "unit circle" at these three values of </>, so there are 
no two-body collisions in this type of orbit. 

Yet, this suggests that there might be other, per- 
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haps multiply periodic solutions with "trajectories" in 
the (r, 4>) plane that touch the r—1 circle at the hyper- 
angle values other than <f> = ±§ and/or approach the 
unit circle (the equator of the shape hemisphere) even 
closer to the two-body collision points. The latter fact 
means that the corresponding trajectories in real space 
are "narrower" than the figure-eight one. 

We have studied this region more closely and found 
several new periodic solutions with lesser symmetry than 
the figure-eight one that pass through the "infinitesimal" 
neighborhood of the initial state, but only in the Y-string 
and the A-string potentials (i.e. not in Newton's gravity, 
as yet). We display two interesting new orbits below. 
The initial conditions are given in Table |TJ 

TABLE I: (Color online) The initial conditions for the solu- 
tions shown in this paper: d is the value of the initial distance 
between the outer left, or right particle and the middle one, 
and the velocities (both are in dimensionless units where the 
masses and the coupling constants have been set equal to 
unity); the angle 9 is in radians. 



name 


d 


V 


#(rad) potential 


fig.8 


6 


1.37 


1.205 Y-string 


type I 


6 


1.32 


1.437 Y-string 


type II 6 


4.53 


1.40 Y-string 


fig-8 


1 


0.6355 


0.5736 Newton 


fig.8 


1 


0.536 


1.49287 A-string 



A. Type I ("linear-in-0") reduced symmetry 
solution 

First note that the real-space trajectory of ("right- 
hand-side") particle number 2 in Fig. [T3]is different from 
the one of the ("central") particle number 1 in Fig. Q~2] 
thus making it clear that this is a periodic, but not a 
choreographic motion. In other words, this solution 
is symmetric "merely" under the two-body permutation 
group S2, rather than under the three-body permutation 
group S3. Due to the reduced symmetry, one particle ex- 
ecutes an "independent" motion, whereas the other two 
move on orbits that are mirror images of each other, very 
much like those in Ref. 18]. This means that this new 
solution is probably a bifurcation of the figure-eight orbit 
as a function of particle masses related to those found in 
the Ref. Q, i.e. as a function of explicit S3 permutation 
symmetry breaking. 

The figure-eight orbit touches the unit circle at three 
points of the equilateral triangle defined by <fi = ± ^tt, it, 
see Fig. El whereas this new solution touches it at only 
one vertex of this equilateral triangle viz. (f> = ^ir, and 
"cuts corners" at the other two, only to touch the unit 
circle at four other values of <fi that are different from the 
two-body collisions points 0, ±§7r. This is still a periodic 
solution with a period of 87r, i.e. it takes four cycles of the 
hyper-angle (f> to complete one period, but with several 
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-8-6 -4 -2 2 4 6 8 

FIG. 12: (Color online) Real space trajectory of particle num- 
ber 1 of the type-I new solution that passes through the figure- 
eight initial configuration. The trajectory of particle number 
3 is a reflection about the line dividing this trajectory verti- 
cally. 
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-1.5 1 1 1 1 1 1 1 1 1 

-8-6 -4 -2 2 4 6 8 

FIG. 13: (Color online) Real space trajectory of particle num- 
ber 2 of the type-I new solution that passes through the figure- 
eight initial configuration. 



different hyper-angular frequencies, instead of the single 
basic frequency 2><f>. This fact may not be immediately 
visible to the naked eye, as these frequencies are close to 
3<p, but shows up as "beats" in the time dependence of 
the amplitudes. 

The hyper-angle <j) in this solution (still) grows (or de- 
scends) indefinitely, so this solution also corresponds to 
a kind of rotation of the pendulum, but with a chang- 
ing angular velocity, see Fig. 1141 The time derivatives 
show the beats more clearly, see Fig. [T5] The temporal 
variation of the hyper-radius R(t) is shown in Fig. [TBI 
and that of the hyper-angular momentum G(t) is shown 
in Fig. [T7] Note the beats in the time evolution of R(t) 
and G(t), as advertised. 
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FIG. 14: (Color online) The time dependence of the hyper- 
angular radius r red (solid), and the hyper-angles a — 
sin -1 r blue (long dashed) and (f> gray (short dashed) of the 
type-I new solution in the Y-string potential that passes 
through the figure-eight initial configuration. The legenda 
are explicitly shown below the figure. Note that cj> moving 
from to 2n corresponds to two segments between vertical 
lines/discontinuities due to the numerical evaluation of in- 
verse trigonometric functions. Note that one complete period 
of the motion corresponds to eight such segments, i.e. to <f> 
changing from to 87r, or to four complete revolutions around 
the (r, <j>) circle. 
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FIG. 15: (Color online) The first derivatives of the time de- 
pendence of the hyper-angular radius f red (solid), and the 
hyper-angles a = 4r sin -1 r blue (long dashed) and <j> gray 
(short dashed) of the type-I new solution in the Y-string po- 
tential that passes through the figure-eight initial configura- 
tion. 



B. Type II ("oscillating 0") reduced symmetry 
solution 

First note that the real-space trajectory of particle 
number 2 in this solution, Fig. [18j is different from the 
one of particle number 1, thus making it clear that this 
is also a periodic, but not a choreographic motion. The 
trajectory of particle number 3 (blue on-line) is a reflec- 
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FIG. 16: (Color online) The time dependence of the hyper- 
radius R of the type-I new solution in the Y-string potential 
that passes through the figure-eight initial configuration. 
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FIG. 17: (Color online) The time dependence of the hyper- 
angular momentum G of the type-I new solution in the Y- 
string potential that passes through the figure-eight initial 
configuration. Note the cut off peaks of the sines ( "flat tops" ) . 

tion of trajectory of particle number 1 about the origin. 
At first, the aforementioned action-angle variables do not 
seem appropriate for this new periodic orbit, indeed the 
(formerly) linear increase/decrease of the hyper-angle <p 
is now subject to substantial modifications: after ini- 
tial rapid (hyper-) rotation in the clock-wise direction 
starting from <f> = 0, it slows down and stops around 
min ~ — 0.7647T, Fig. [191 then changes the direction of 
motion and swings back yet again only to stop, this time 
around </> max ~ 1.4317T, then repeating this cycle ad in- 
finitum. Note that the maximal difference (twice the am- 
plitude) of cj) is numerically close to being a simple frac- 
tion of 7r, i.e. A(f) — </> max — </> m i n = 13.0001^, whereas 

the average value <j) = |(<Anax + 0min) is numerically close 
to j7r. We suspect that ^ir and \tt are the exact val- 
ues and that the deviations from our numerical values 
are due to rounding-off errors. This resembles the os- 
cillations of a variable-length pendulum. Indeed, FigfTOl 
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-8-6 -4-2 2 4 6 8 

FIG. 18: (Color online) Real space trajectories of particle 
number 1 red (solid), particle number 2 green (light gray 
dashed) and particle number 3 blue (dark gray dashed) line in 
the type-II new solution that passes through the figure-eight 
initial configuration. The legenda are explicitly shown in the 
lower right-hand side corner of the figure. 

60 I 1 1 1 1 1 1 1 1 1 
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R 10o 

FIG. 19: (Color online) The time dependence of the hyper- 
radius R red (solid), and 10x0 gray (short dashed) in the 
type-II "oscillating" solution in the Y-string potential. The 
legenda are explicitly shown below the figure. 

shows that the hyper-radius R(t) is oscillating with the 
same frequency and phase as the hyper-angle <j)(t), thus 
extending the analogy with the variable-length pendulum 
model. 



V. CONCLUSIONS 

We have studied the figure-eight motion in three differ- 
ent three-body potentials in terms of permutation sym- 
metric variables (R,r = sin a), and the braiding hyper- 
angle 4>. The existence of this orbit depends on the peri- 
odic dependence of the potential on the braiding hyper- 
angle 4> that "guides" the figure-eight orbit (s) around the 



two-body collision points. 

The figure-eight orbits in the triangle shape space are 
generally close to their iso-potential lines, though formal 
arguments show that they cannot be exactly identical 
4] . Thus the exact analytic solutions ought to be sought 
among (small) oscillations about the iso-potential lines, 
with basic frequency 30. 

The Hamiltonian of three identical particles in a per- 
mutation symmetric potential with vanishing total an- 
gular momentum has certain similarities with that of 
a spherical pendulum in inhomogeneous azimuthally- 
periodic potentials, which, in turn, suggests existence of 
other types of solutions. 

We have found two new periodic solutions in the Y- 
string potential that pass through the figure-eight initial 
state, but do not share its symmetry. One of these so- 
lutions (type I) has a monotonically rising/ descending 
hyper-angle 4>, just like the figure-eight orbit, but a dif- 
ferent pattern of syzygies, whereas the second (type II) 
new solution's </> is oscillating about its average value of 
^, with the hyper-radius R following suit. 

All of these orbits are clearly characterized by their 
{R, r, (f>) behaviors that display certain similarities, de- 
spite their independent and seemingly random form of 
trajectories in the configuration space. Thus we believe 
this to be a good set of variables to mathematically sim- 
plify and describe all periodic orbits of three identical 
bodies. 

A few words about the history of this subject and our 
approach to it might be in order now. In our study Ref. 
[1J| of the so-called Y-junction and the A-string poten- 
tials we found an integral of three-body motion, when 
the three-body potential depends on only two, rather 
than three, independent three-body variables, viz. the 
hyper-radius (or the moment of inertia divided by the 
quark/particle mass) and the area of the triangle defined 
by the three bodies [23j . As the static three-body poten- 
tial may depend on (at most) three independent scalar 
variables, our observation naturally begged the question: 
what is the third independent three-body variable in this 
set? [HJ It was an attempt to answer this question that 
brought us to the present permutation symmetric vari- 
ables. We are not the first ones to use them, however: 
Chenciner and Montgomery have used these variables 
(these authors call 9 "our" variable 4>) to parametrize 
the triangle shape space in Ref. 0]. According to Ref. 
0,H. Hopf was the first one to introduce these variables, 
Ref. [H. 
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